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Abstract. We investigate the behaviour of cosmological baryons at low redshifts z £ 5 after reionization through 
analytic means. In particular, we study the density-temperature phase-diagram which describes the history of 
the gas. We show how the location of the matter in this (p, T) diagram expresses the various constraints implied 
by usual hierarchical scenarios. This yields robust model-independent results which agree with numerical simula¬ 
tions. The IGM is seen to be formed via two phases: a “cool” photo-ionized component and a “warm” component 
governed by shock-heating. We also briefly describe how the remainder of the matter is distributed over galaxies, 
groups and clusters. We recover the fraction of matter and the spatial clustering computed by numerical sim¬ 
ulations. We also check that the soft X-ray background due to the “warm” IGM component is consistent with 
observations. We find in the present universe a baryon fraction of 7% in hot gas, 24% in the warm IGM, 38% in 
the cool IGM, 9% within star-like objects and, as a still un-observed component, 22% of dark baryons associated 
with collapsed structures, with a relative uncertainty no larger than 30% on these numbers. 
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1. Introduction 

As is well known, the mass of baryons embedded within 
stars or galactic disks in the current universe is quite 
small since it yields a baryonic parameter fl ga i ~ 


0.002 — 0.006 (e.g., Fukugita et al. (1998)) while stan¬ 
dard nucleosynthesis calculations give fib — 0.045 (e.g., 
Tytler et al. (2000)). Therefore, most of the baryonic mat¬ 
ter should lie in the intergalactic medium. This agrees 
rather well with the fact that at higher redshift one 
observes a large amount of mass in the Lyman-a for¬ 
est which consists of moderate density fluctuations ion¬ 
ized by the background UV flux emitted by distant 
galaxies. Thus, one gets fRya ~ 0.01 — 0.05 at z 


3 (e.g., |Fukugita et al. (1998)[) . However, as noticed in 
Cen fc Ostriker (1999) the mass within the Lyman-a for¬ 
est decreases with time and at z = 0 summing over all 
observed contributions one obtains fib ^ 0.011 which 
falls short of the required total baryonic mass. Hence 
at z = 0 a large part of the baryons must lie in a 
new intergalactic component beyond the usual Lyman- 
a forest clouds. As argued in Cen & Ostriker (1999) and 
Dave et al. (2001)| this could be part of a “warm” phase 
of the intergalactic medium (IGM), with temperatures in 
the range 10 5 < T < 10 7 K. 

The latter conclusion was reached from numerical sim¬ 
ulations. In this article, we reconsider this problem in or¬ 
der to derive the properties of the IGM by analytic means. 


In particular, we wish to investigate whether one can un¬ 
derstand this behaviour in a quantitative manner from 
robust, model-independent, arguments. 

First, in Sect.|| we study the (p, T) phase-diagram 
of cosmological baryons. While the Lyman-a forest 
is described by a well-defined Equation of State the 
“warm” IGM component shows a broad scatter (e.g., 
Dave et al. (1999)| , Dave et al. (2001)) since its temper¬ 
ature depends through shock-heating on the neighbour¬ 
ing gravitational potential which is a stochastic field. 
Nevertheless, we show that it is constrained to lie in a 
well-defined domain in the (p, T) plane, and determine its 
average location in this plane, which may be considered 
as the “Equation of State” of the “warm” IGM. We also 
give the location in this diagram of galaxies, groups and 
clusters. 

Next, in Sect.|| we use our results to compute the red- 
shift evolution of the fraction of matter enclosed within 
the different phases. Then, in Sect.|| we estimate the two- 
point correlation function and the clumping factor of the 
“warm” IGM. Finally, in Sect.|5| we check that the X- 
ray background emitted by the “warm” component agrees 
with observations. 


2. The phase-diagram of cosmological baryons 

In this section we investigate the evolution of the physi¬ 
cal properties of the baryons at moderate redshifts z £ 5, 
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as gravitational clustering builds up. To do so, we study 
the distribution of matter in a density-temperature plane 
and determine the (p, T) phase-diagram. Indeed, as we de¬ 
scribe below baryons are located within specific regions in 
this (p, T ) plane which reflect the various processes (cool¬ 
ing, radiative heating, shock-heating, gravitational clus¬ 
tering,...) which take place in the universe. We display 
our results in Fig.|l]. Note that we use the overdensity 
l + <5 = p/p rather than the density p in the figures, where 
p is the mean density of the universe. 

In this article, all our numerical results are obtained 
with the following cosmological parameters. We consider 
a low-density flat universe with = 0.4 and = 0.6. 
The baryonic density parameter is Sib = 0.0473 and the 
Hubble constant is Hq = 65 km/s/Mpc. The ACDM 
power-spectrum of the linear density fluctuations is nor¬ 
malized by as = 0 . 8 . 


2.1. Temperature-scale relations 

We consider two processes for the possible heating of the 
gas. First, the gas temperature may be governed by a lo¬ 
cal shock-heating due to gravitational processes. Thus, as 
large-scale structures build up we expect tidal effects and 
virialization processes to heat the gas within and in the 
neighbourhood of collapsed halos and filaments up to the 
virial temperature of the underlying gravitational poten¬ 
tial well. The latter scales as $ ~ QpR 2 , which defines the 
characteristic scale R, hence we write: 


local heating : 


kT = —pm p GpR 2 , 


(1) 


where pm p is the mean molecular weight of the gas and 
Trip the proton mass. 

Second, the gas temperature may be set by an external 
heating source, like a UV background radiation. Then, 
the temperature is no longer given by eq. ( 0 ) since it is not 
fixed by gravitational energy. However, the pressure of the 
gas erases the baryonic density fluctuations over a scale R 
given by: 


external heating : R = tnC s 



( 2 ) 


2.2. Exclusion constraints 

First, we note that the distribution of matter in the (p, T ) 
plane is strongly constrained by the existence of some ex¬ 
clusion regions where no particles can be found. The ad¬ 
vantage of such conditions is that they are very robust in 
the sense they do not depend on the detailed history of 
the baryons. They apply to any hierarchical scenario of 
structure formation and provide the first guideline for the 
properties of baryons. We discuss these various constraints 
in the sections below. For sake of clarity, we turn within 
the diagrams shown in Fig.^ in counter-clockwise order. 


2.2.1. Cooling constraint 


The dashed-line on the far-right with an inclined “V- 
shape” is the familiar cooling constraint (see also 


Rees fe Ostriker (1977% pilk (1977)|) . It expresses the fact 
that high-density ionized hot gas (T > 10 4 K) cools very 
rapidly. Thus, baryons cannot remain for long in regions 
located to the right of this dashed-line. More precisely, 
we define the “effective” cooling time t c ,e f set by atomic 
physics as: 


cooling+UV 


^cool ^heat 


e 

^c,ef 


= - —, (3) 


where e = 3/2 ribkT is the specific energy and rib is 
the baryon number density. The cooling time t C ooi de¬ 
scribes collisional excitation, collisional ionization, recom¬ 
bination, molecular hydrogen cooli ng, bremsstrahlung an d 
Compton cooling or heating (e.g., Anninos et al. (1997) ). 
The heating time fheat is the radiative heating time due 
to a UV background radiation of the form J „ oc u~ l . 
As usual, we normalize the UV flux by its value J 21 
at the HI ionization threshold (912A) in units of 1CU 21 
erg s _1 Hz -1 cm -2 sr _1 . The values we choose for J 21 
are given in Tab.l as a function of redshift. They are 
consistent with observations (e.g., piallongo et al. (1996) 


|Cooke et al. (1997)|, Vogel et al. (1995)) and are the same 


as the ones used for modeling the evolution of Lyman- 
a clouds Valageas et al. (1999) . Then, the cooling curve 
shown in Fig.|lJ is given by the condition: 


cooling : t c , e f = in- 


(4) 


where C s is the sound velocity and 7 = 5/3. This is simply 
the distance over which sound waves can travel in a Hubble 
time and build pressure equilibrium. Therefore, the gas 
sees the dark matter density field smoothed over this scale 

R. 

There is a characteristic density p c at which shock 
heating starts to play a role, that is where eq.(jl) is to 
be used in place of eq.(|J). It is expected to correspond 
to a density contrast of a few units at least. The detailed 
discussion of the relation bet ween the two regimes around 
p ~ p c may be found in Sect.Td below. 


Note that in eq.([|) we neglected shock-heating due to 
gravitational processes (e.g., virialization) and possible 
adiabatic cooling due to the expansion of the fluid ele¬ 
ments. However, these processes are governed by the time 
tn which measures the expansion rate of the universe and 
the time-scale over which new non-linear structures form. 
Therefore, they do not modify the form of eq.([|): they 
may only change the ratio between both sides of eq.fji]) 
(which we set equal to one) by a factor of order unity. 
Moreover, we have f Cje f — i coo i as radiative heating does 
play a strong role here. Hence this constraint (0) is actu¬ 
ally quite robust. Then, the regions located to the right of 
the cooling curve have a cooling time which is less than 
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Table 1 . Redshift evolution of the background UV flux 
J 2 i(z) used in this article. 


z=0 1 2 3 4 5 


J 21 = 0.05 0.5 0.8 0.4 0.2 0.1 


the Hubble time. As a consequence, no particle can re¬ 
main in this part of the (p, T ) plane for a Hubble time. 
This process is actually at the origin of galaxy formation 
as this gas cools and falls to the center of the gra vita¬ 
tional well to form a disk and stars (see Sect.2/7 and 
Valageas fe Schaeffer (1999) ). This means that the matter 
which entered this region at some time should accumulate 
on the lower branch of the cooling curve (at T ~ 10 4 
K) under the form of high-density cool gas which evolves 
on a Hubble time (but some of this gas also turns into 
stars). Note that this accumulation of baryons along the 
lower branch is clearly seen in the numerical simulations 
displayed in Fig.11 in Dave et al. (1999) and in Fig.9 in 


Springel fe Hernquist (2001) 


2.2.2. High-density fluctuations 


A second constraint is set by the properties of the dark 
matter density field itself. Indeed, at a given scale R be¬ 
yond some threshold p+(R) large densities have a negli¬ 
gible probability and only involve a very small amount of 
matter. In the linear regime for instance, at large scales 
R } the probability distribution function (pdf) V(Sr) of 
the density contrast Sr over a spherical cell of radius R 
is a Gaussian of variance a(R). Here we note as usual 
cr(R) the rms linear density fluctuation and we assume 
Gaussian initial density fluctuations. In the non-linear 
regime, for highly non-Gaussian probabilities, the ana¬ 
logue of this upper threshold can also be calculated. It 
is given in App.[A] where we detail our model for the pdf 
P(p R ), where pR = 1 + Sr is the overdensity. The high- 
density cutoff defined by eq.(A.E) obeys the asymptotic 
behaviour (A. 10) in the quasi-linear and highly non-linear 
regimes. Note that this asymptotic behaviour actually is 
model-independent. Indeed, in the limit £ -C 1 we re¬ 
cover the Gaussian cutoff implied by the initial conditions 
(where £ = (Sr)). On the other hand, in the limit £ 1 

we recover the fact that most of the matter is enclosed 
within high-density halos of density contrast Sr ~ £ which 
occupy a small fraction of the volume (~ l/£) of the uni¬ 
verse. Therefore, the constraint given by the high-density 
cutoff of the pdf V(pr) in the (p, T ) plane is quite robust 
and it applies to all hierarchical scenarios of structure for¬ 
mation. 

We assume that the baryon density p& scales as the 
dark matter density through 


Pb = 


Db 


(5) 


This is valid even in the cooling region discussed above in 
Sect. 2.2.1 if we consider also the cold baryons. We then can 
obtain an upper cutoff for the baryon overdensity p+ ( R ) at 
any given scale R. It will turn out to be convenient to ex¬ 
press this cutoff p+(R) in terms of a density-temperature 
relation, in order to draw its consequences for the phase- 
diagram ( p , T) of the IGM. To do so we simply need the 
relation T(p, R). 

Substituting eq. ( 0 ) into the relation p+(R) derived in 
eq.(A.S) we obtain a curve p+,i oc (T). The subscript “loc” 
in the overdensity cutoff p+ i oc refers to the fact that 
this is a “local” heating process. It is due to the gravita¬ 
tional interaction with neighbouring structures. This high- 
density cutoff ioc(T) corresponds to the dot-dashed 
curve shown in the diagrams in Fig.|l] which runs from 
S ~ 5 up to S ~ 10 4 and which crosses the cooling curve. 

Subst itut ing eq.(|j) into the relation p+(R) obtained 
from eq.( |A. 8 | ) we get a curve p+, ext (T). This corresponds 
in Fig.p] to the short branch which runs upward from 
the curve p +i i oc (T) at S ~ 10. Indeed, we note that for 
<5 10 the gas is located close to non-linear structures 

so that local shock-heating must be taken into account. 
However, for large densities we have T\ oc > T ext since 
eq.(j|) yields T ext ~ G~pR 2 . In this case external heating 
plays no role. This is why we only plot the curve p+, ex t(T) 
up to the characteristic density contrast where Xi oc = Text. 
At higher densities the threshold p+ iex t(T) becomes irrel¬ 
evant. In fact, we see in Fig.jl] that the curve p+, ex t(T) 
due to external heating plays no role since at moderate 
densities p ~ p it is repelled to quasi-linear scales (see 
Sect . [2. 2. 3[ ). 

Then, the region to the upper-right of the curves 
p+,ioc (T) and p+, ext (T) in the (p, T) plane shown in Fig]l] 
corresponds to rare high-density fluctuations which are lo¬ 
cated in the tail of the pdf V(pr). Therefore, there should 
be very few particles beyond these lines. Thus, this defines 
a second exclusion region. 

Note that at z = 0 this constraint yields an upper 
bound T+ ~ 10 7 3 K ~ 1.7 keV. Of course, there exist 
some halos with a larger temperature: massive X-ray clus¬ 
ters. However, these are rare objects which only contain a 
small fraction of the baryonic matter content of the uni¬ 
verse (typically 10 %) and they indeed correspond to the 
high-mass tail of the mass function. This is obviously con¬ 
sistent with the description of the baryonic matter which 
is worked out in this article, see Sect.2.7. 


2.2.3. Transition to non-linearity 

A characteristic scale which enters the problem we investi¬ 
gate here is set by the transition to the non-linear regime. 
At a given epoch, we define this scale Rq(z) by the rela¬ 
tion: 

<t(Rq,z) = 1. ( 6 ) 

Thus, scales larger than Rq are still within the linear 
regime. Then, there can be no shock-heating due to grav¬ 
itational clustering on these large scales which have not 
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turned non-linear yet. Using eq. ( 0 ) we obtain a character¬ 
istic temperature T n \(p): 

87T 

transition to non-linearity : kT n \ = — pm p QpR^ } (7) 

which describes the transition to non-linear scales (which 
lie at T -C T n i(p)). This is shown by the straight dashed- 
line (T n i oc p) plotted in Fig.p]. Besides, we note that 
known sources of external heating (e.g., the UV back¬ 
ground radiation) cannot heat the IGM up to such high 
temperatures. Therefore, no fluid element can be located 
in the region to the upper-left of the curve T n \(p) since 
no physical process which is active on these large scales 
can heat the gas to these high temperatures. This yields 
a third exclusion region. 


2.2.4. Low-density fluctuations 

Finally, a fourth constraint on the distribution of mat¬ 
ter is given by the low-density cutoff of the pdf V(pn) 
of the dark matter density field. This i s the analog of 
the high-density cutoff discussed in Sect . [2. 2.2 . The low - 
density cutoff p~(R) is derived in App.A, from eq.(A.8). 
It obeys the asymptotic behaviours ( A.11 ). In the quasi- 
linear regime we again recover the usual Gaussian cut¬ 
off S- ~ —a. In the highly non-linear regime the small 

overdensity p- ~ £ <C 1 expresses the formation 
of extreme underdensities on small scales. Contrary to 
the high-density cutoff it is somewhat model-dependent 
through the exponent k but this has no strong effect on 
the (p, T) ph ase-diagram (we typically have n ~ 0.8). As 
in Sect . [2.2.2 we need a relation T(R) in order to derive 
a condition of the form p-(T). We again consider both 
cases of local and external heating, described by eq.(l) 
and eq.(^j). This yields the curves p_^ oc (T) and p_ jext (T) 
shown by the two steep parallel dashed-lines in Fig.|l], at 
Pr ~ 10 -1 . The curve associated with external heating is 
the left one (i.e. lower densities or higher temperature), 
as can be seen from eq. (ill) and eq.(|2]). Thus, this defines 
a fourth exclusion region to the left of these curves. 


2.3. Equation of State of the Lyman-a forest 

Thus, so far we have obtained constraints on the distribu¬ 
tion of matter in the (p, T) plane by drawing four exclu¬ 
sion regions. This already gives quite useful information 
about the properties of the IGM which are very robust. 
Now, we investigate a different point, seeking the location 
of the gas in the (p, T) plane. This amounts to deriving 
an Equation of State for this component. 

As shown in [I lui fe Gnedin (1997) the low-density 
photo-ionized IGM exhibits such an Equation of State as 
the gas follows a specific relati on T a (p) with a rather small 
scatter. This was derived in Hui fe Gnedin (1997) from 
the Zel’dovich approximation ( Zel’dovich (1970) ) which 
applies up to the moderately non-linear regime (£ ^ 1). 
Here we reconsider this problem and we show that this 
Equation of State is rather robust with respect to the past 


history of the gas and applies independently of the validity 
of the Zel’dovich approximation. First, we assume photo¬ 
ionization equilibrium (we restrict ourselves to z ^ 5 af¬ 
ter reionization) and we only take into account Hydrogen. 
Therefore, the ionization equilibrium reads: 

T rini = ol(T) nmin e , (8) 

where T is the photo-ionization rate and a(T) is the re¬ 
combination rate. They are given by: 

/ i 

47tJ„o-hi^ = 3.08 x 10~ 12 J 2 i(z) s" 1 (9) 

and: 

/ t* \ — i) 

a(T) = aj ^ — J , ( 10 ) 

where we defined: 


v = 1.7, aj = 1.23 x 10 -13 cm 3 s _1 , Tj = 5.8 x 10 4 K.(11) 

The temperature Tj we introduced in eq. © is the char¬ 
acteristic temperature reached by the gas through the 
heating due to the UV background radiation flux. It is 
given by: 


J nvru^c/HI hv 


( 12 ) 


where hvni = 13.6 eV is the Hydrogen ionization thresh¬ 
old. Note that this temperature Tj does not depend on the 
amplitude J 2 i of the UV background. Moreover, it is fixed 
by atomic physics, independently of cosmological param¬ 
eters. Next, the temperature T of a given fluid element 
evolves as: 


1 dT _ 2 1 dp 1 
T dt 3 p dt + fheat 


(13) 


where d/dt is the Lagrangian time derivative. The heating 
time theat is given by: 


3/2nbkT 

the&t 


J dirJ v crmnu.\{hv - hv ffl ) 


dv 

hv 


kTjmul\ 


(14) 


where we take rib = 2 pb/m p since we approximate the gas 
as fully ionized Hydrogen. The evolution eq. © is the 
same as eq.(||), except that we neglect cooling (which is 
justified here since we consider here moderate densities 
and temperatures) and gravitational shock-heating which 
is irrelevant. On the other hand, the term dp/dt represents 
the pressure work, which takes into account the expansion 
of the fluid element. Substituting eq. & into eq.(|T3|) we 
obtain: 


1 dT _ 2 1 dp 
T dt 3 p dt 

where we define: 
pj = 3 


J __P_ 

tio Pj 


-V 

Tj) 


Ub ajt io 


~ 1.1 x 10 28 g cm 3 


(15) 


(16) 
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and: 


f 10 = 10 10 years. 


(17) 


In eq.(^5|) we used the fact that Hydrogen is almost en¬ 
tirely ionized at the low redshifts which we consider here 
(z < 5, after reionization), as shown by the Gunn-Peterson 
test. Next, it is convenient to introduce the adiabat I\ de¬ 
fined as the exponential of the specific entropy s. More 
precisely, we define K and s by: 


A ' = li(£f /S ' * sl " A 

Then, eq. © may be written: 


K 


,-id K 
d t 


1 

Ho \pj 


1-2i//3 


(18) 


(19) 


Here we note that 1 — 2v/3 ~ —0.13 is a small number. 
Therefore, the specific entropy of the fluid shows a weak 
dependence on the evolution of its density p. Hence we 
approximate the solution of eq. @ b y : 

/ \ 1 — 2i//3 , 

K V — V (20) 

\PJJ Ho 


where p is the density at the time t we consider. In terms 
of the temperature-density relation, this yields: 


T a (p) =Tj[v 


*_P_ 
Ho Pj 


l/v 


( 21 ) 


where we used eq. ( |l8| ) . Thus, the expression gives the 
first-order term for the expansion of In T in terms of the 
small parameter (1 — 2i//3)(dlnp/dlnf). From eq. ( 0 ) we 
note that at z ^ 5 we have p ~ pj and tn ~ Ho hence 
the temperature of the photo-ionized IGM will be of order 
Tj ~ 10 4 K. 

Note that eq.(|2i]) is independent of the normalization 
of the UV flux J„. Indeed, the efficiency of radiative heat¬ 
ing is proportional to J„nm but the density of neutral 
Hydrogen scales as 1/J„ (at ionization equilibrium for al¬ 
most fully ionized gas) so that J v cancels out. Therefore, 
the result (^l|) is quite robust since it does not depend 
on the value of the UV flux. In particular, eq. (fill) still 
holds even if the UV background is inhomogeneous: we 
only need to assume local ionization equilibrium. This ex¬ 
plains why there is only a very small scatter around the 
Equation of State (|2l]) since the actual physical conditions 
within these small clouds are almost independent of the 
actual history of each fluid element (i.e. the evolution of 
its density and local UV flux). This indeed agrees with the 
results of numerical simulations (e.g., pave et al. (1999)| ). 

We display in Fig.[l] the Equation of State (|2l|) as 
the solid line which runs through T~10 4 Kat(i = 0. 
Beyond the overdensity p c defined in eq.(|2|) below, we 
plot this line as a dashe d-line , until it enters the cooling 
regi on des cribed in Sect. 2.2.1 . Indeed, as we explained in 
Sect. 2.2.2 at large densities the “virial temperature” Ij oc 
becomes larger than the temperature T ext due to some 


external energy source (here photo-ionization heating by 
the UV background radiation). This means that for these 
regions, which have already reached the non-linear regime 
as 5 > S c ~ 5 (see eq. ©), shock-heating due to the 
gravitational dynamics can no longer be neglected and 
it actually becomes dominant. Therefore, for high densi¬ 
ties with p > p c the gas should no longer fall onto the 
curve ©■ Nevertheless, since shock-heating can only in¬ 
crease the temperature of the gas the relation (^l]) now 
provides a lower bound to the temperature T. Hence, the 
region below the curve (21) is excluded in the (p,T) di¬ 
agram. This holds until we enter the cooling region dis¬ 
cussed in Sec t. 2.2.1 . On the low-density side, as described 
in Sect . [2. 2 .4 we are constrained by the low-density cutoff 

p~ : e X t(T). 


Therefore, we predict that we should have two phases 
for the IGM. A first “cool” phase is described by 
the Equation of State © with intermediate densities 
P-,ex t < P < Pc- It is photo-ionized gas heated up to 
T ~ 10 4 K by the background UV flux. This corre¬ 
sponds to the moderate density fluctuations which form 
the Lyman-cc forest. A second “warm” phase is made 
of higher-density regions which have already experienced 
some shock-heating due to the building of gravitational 
structures but which have not entered the cooling region 
yet. These particles should be located in the (p,T) plane 
above the curve © and within the constraints described 
in the previous sections. 


2.4. Warm IGM 


Thus, the previous discussion shows that a “warm” com¬ 
ponent of the IGM should no longer follow the equation 
of state (|2l]). It however remain s co nstrained by the ex¬ 
clusion regions obtained in Sect. 222, This leaves a closed 
allowed region in the (p, T) plane. Since the temperature 
of the gas should be governed by the shock-heating due to 
the gravitational dynamics we can expect a rather large 
scatter. Indeed, the history of these particles depends on 
the properties of their neighbouring dark matter density 
field. This introduces a stochastic component which gives 
rise to a wide variety of possible histories so that there 
is no longer a unique temperature-density relation. This 
agrees well with Fig.11 in Dave et ah (1999) which shows 
indeed a broad cloud of points in the (p, T) plane within 
the allowed region we obtained in the previous sections. 
Note also that at z = 0 some of the gas which has been 
shock-heated to T > 10 4 K actually lies at low density 
contrasts 6 < 0 in the numerical simulation. This is con¬ 
sistent with our results shown in Fig|I|. This simply means 
that on small non-linear scales even regions with S ^ 0 
have experienced shell-crossing and gravitational shock¬ 
heating (this also implies that they are not just described 
by a collapse dynamics). 

However, the distribution of matter obtained in 


Dave et al. (1999) from numerical simulations does not fill 
entirely this allowed region (note indeed that this is not 
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0 3 4 


log(l+<5) 


Fig. 1 . The phase-diagram of the IGM from 2 = 3 down to z = 0. The straight solid line T a shows the Equation of 
State of the “cool” IGM (Lyman-a forest). The curved solid line T g h shows, as a mean trend, the “Equation of State” 
of the “warm” IGM which is shock-heated through the building of non-linear gravitational structures. The dashed 
curves draw exclusion regions around the allowed domain for this “warm” IGM. In counter-clockwise order, starting 
from the right side, they correspond to 1) a fast-cooling region where the gas cannot remain over a Hubble time, 2) a 
high-density and high-temperature domain within the exponential tail of the pdf V{5r) which only contains very rare 
massive halos, 3) large linear scales where gravitational shock-heating has not appeared yet, 4) a low-density region 
within the tail of the pdf V(5r) which is associated with very rare voids and 5) a lower-bound for the temperature T a 
set by radiative heating from the UV background. The vertical dotted line is the density threshold A V j r of just-virialized 
halos. The points show the results of numerical simulations from Dave et al. (1999) (Fig.11). 


implied by our previous considerations) and it is not uni¬ 
form. Hence it would be convenient to derive a curve in the 
(p, T) plane which would describe the “mean” behaviour 
of this warm phase of the IGM. To this purpose, we first 
compute the gas temperature associated with a given halo 
of mass M which forms at redshift z. As noticed above, 
this does not apply to low-density regions with 6 ^ 0 but 
this model should provide a useful estimate for the higher 
density regions which actually contain most of the mass. 
Gravitational heating is effective in the collapsing phase 
of the halo which may be set at 5 > 5 C . One expects S c to 
be of the order of a few units. We set it to 

1 + S c = 6 and p c = (1 + <5 c )p, (22) 

which is close to the value given by the spherical collapse 
dynamics at turn-around. We consider this not to be an 


accidental coincidence. Indeed, the kinetic energy of the 
underlying overdense large-scale structure is not thermal- 
ized during its expansion phase and it only provides an 
energy source for the heating of the gas within the collapse 
stage. Even for an expanding cloud which is not perfectly 
spherical, such a phase of maximum expansion is reached, 
where the kinetic energy in the collective motion is small. 
Afterwards, the cloud enters a collapse phase where shocks 
develop. As is usual for shocks, the latter transform col¬ 
lective large-scale kinetic energy into small-scale thermal 
energy. The collective kinetic energy per particle v 2 in the 
collapse phase (i.e. the squared velocity of the inflow) is of 
order v 2 ~ [QM/R — i>q] where the constant term Vq ex¬ 
presses the fact that at turn-around v 2 = 0. Then, this ki¬ 
netic energy is locally transformed into heat in the shocked 
regions. The temperature reached there T g h is of the order 
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of this kinetic energy since it is set by the inflow velocity v 
through the matching conditions at the interface between 
the inner warm shocked region and the outer collapsing 
region. We thus obtain a relation between the small-scale 
temperature (which applies to the small shocked regions) 
and the large-scale kinetic energy per particle. This yields 
(for a large-scale structure of mass M ): 


T gh 


T c 



(23) 


with 

Tc = |(f) 1/ %m p »M 2/ y; / 3. (24, 

Of course, the characteristic temperature T c depends on 
the mass M of the collapsing cloud. The subscript “gh” 
refers to the fact that the temperature T g h is governed by 
“gravitational heating”. This provides an estimate T g h(p) 
of the location of the warm spots in the (p, T) plane, for 
fluid elements parameterized by the mass M of the un¬ 
derlying larger-scale structure which is collapsing. Here 
we assumed that the density within the shocked regions 
scales as the overall density of the collapsing structure 
(note that a strong shock only increases the density by a 
factor 4 at most). 

Now, seeking an average location of these warm spots, 
we estimate the typical mass M 0 (z) which collapses at 
each redshift z by: 


o-(M 0 , z) = 1. (25) 

This also yields the scale Rq(z) defined in eq.(p|), with 
M 0 = 47r/3pf?Q. Among all the T g h(p) relations (^3[) which 
have been seen to depend on the halo mass M, we can then 
find a “mean” locus by taking a characteristic temperature 
T c , eq.((2j), calculated with the mass M 0 . 

The mass Mq and the scale Ro characterize the non¬ 
linear structures which just turned non-linear. For in¬ 
stance, at z = 0 we have Rq ~ 9 Mpc. At virial- 
ization when the density contrast reaches the threshold 
A v ; r ~ 200 this yields a radius R vir ~ 1.3 Mpc as for 
large clusters. On the other hand, the shocked regions, at 
temperature T g h, correspond to much smaller scales. Their 
size R can be estimated by using approximate hydrostatic 
equilibrium within the filaments along the transverse di¬ 
rection^]. As a consequence, the eq.(^) still holds, up to 
factors of order unity, if R is taken to be the thickness of 
the filament. This is a natural result since the inflow and 
the shock actually are at the origin of the filament. Also, 

1 Indeed, the pressure beyond the shock is of order P ~ pT g h 
while the velocity drops. The Euler equation then yields P/p ~ 
<1> where <t> is the local gravitational potential well associated 
with the filament. Note that we do not need exact hydrostatic 
equilibrium within this shocked region: e.g. there may be a non¬ 
zero flow along the axis of the filament. We only use the fact 
that the kinetic energy associated to the transverse collective 
velocity is small as compared to the local pressure because the 
shock converts the infall velocity into thermal energy. 


the thickness R of the filaments turns out, as a rule, to be 
much smaller than Rq. For instance, at z = 0 we find it 
is of order of a few hundred kpc, see Sect.|] below. Note 
that the same processes are at work for the dark matter 
density field, except that shocked regions now correspond 
to areas where shell-crossing governs the dynamics and 
builds a large velocity dispersion. Therefore, we can again 
assume that the gas follows the dark matter density field 
and the gravitational potential is dominated by the dark 
matter. Our findings agree with the results of numerical 
simulations which show that the shocked regions are asso¬ 
ciated with the filaments which appear in the dark matter 
density field. Thus, as sho wn by the simulation map in 
Fig.3 in Dave et al. (2001) , the “warm” IGM is an intri¬ 
cate network of filaments which extends over a few Mpc 
(the scale Rq) but the thickness of the filaments (the scale 
R) is much smaller than this global scale. In Sect.[| we will 
see in addition that the scale R of the warm regions shows 
up as a small-scale cutoff for the “warm” IGM two-point 
correlation function. We obtain a length of a few hundred 
kpc at z = 0 which is much smaller than the few Mpc 
which characterize the underlying global structure which 
is collapsing. This also agrees with the numerical simula¬ 
tions of Dave et al. (2001)| . 

Therefore, the “warm” IGM component is described 
by the curve (^3|) in the (p,T) plane, together with eq. (JTJ) 
which yields the local size R of the clouds (or the thickness 
of the filaments). Our model relies on three major points: 

1) there is a critical overdensity associated to the “warm” 
IGM which reflects the turn-around of patches just going 
non-linear, 2) shock-heating locally transforms into heat 
the collective kinetic energy and 3) the scale of the shocked 
regions is given by local approximate hydrostatic equilib¬ 
rium at the latter temperature. This yields 1) the locus in 
parameter space (p, T) where the “warm” IGM appears, 

2) the size of the non-linear network of shocked regions 
(filament network) which is the scale Rq just turning non¬ 
linear and 3) the thickness R of the filaments which ap¬ 
pears to be much smaller than the previous scale. 

We also see that at the beginning of the cloud collapse 
the shocked regions are predicted to have low tempera¬ 
tures correlated with a rather small spatial extension. As 
the overall collapse proceeds their temperature and their 
size increase and the highest temperatures are reached 
when the cloud is being virialized with scales of the order 
of the radius of the halo. This is more or less what could be 
expected to occur. Thus, we model the mean trend of the 
“warm” IGM by eq.(|2^) which appears as a curved solid 
line labeled ‘Tgh” in Fig|j]. It starts from the end-point 
of the “Cool” IGM (at <5 ~ 5) and it goes towards larger 
densities and temperatures until it reaches the cooling re¬ 
gion (at z 3) or the density threshold A v j r (at low z). 
Indeed, beyond this point we consider that we have cooled 
objects (galaxies) or hot virialized halos (clusters) which 
are not part of the “warm” IGM. 

As seen from Fig.0. at low z gravitational heating 
rapidly dominates over the UV heating, for densities close 
to p c and temperatures around T c = T a (p c ). On the other 
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hand, the “warm” IGM clouds are not the same as the 
UV heated objects which form the “cool” component. 
Note that the curve (|2^) used with the mass (^||) only 
reflects the average trend. Indeed, the stochastic char¬ 
acter of the dark matter density held leads to a broad 
variety of masses which are just collapsing, and hence 
of (p,T) trajectories. Moreover, the local properties of 
the shocks also provide for some additional scatter. This 
induces a (rather large) dispersion of the points of the 
“warm” IGM in the (p, T ) plane, which agrees with nu¬ 
merical simulations (e.g., Dave et al. (2001)). However the 
“warm” IGM should remain enclosed within the allowed 
region defined by the constraints discussed in the previ¬ 
ous sections. Indeed, we must point out that the valid ity of 
eqJ3) ensures that the constraints obtained in Sect . |2.2.2) 
- 2.2.4 still hold. 


2.5. The different phases of the IGM 

In the previous sections we have shown that one can distin¬ 
guish two components in the IGM and we have determined 
their location in the (p, T) phase-diagram. Our results are 
displayed in Fig.[lj. 

Firstly, there is a “cool” IGM phase (T ~ 10 3 — 10 4 
K) which corresponds to the Lyman-a forest. These are 
moderate density fluctuations (<5 < 5 C ) of photo-ionized 
gas. They are described by the Equation of State ( |2l| ) 
which arises from the heating of the gas by the UV back¬ 
ground and the cooling due to the expansion (i.e. pressure 
work). Thus, this component lies on a well-defined curve 
in the (p, T) plane. This curve is bounded towards low 
densities by the cutoff of the pdf V(pr ), which expresses 
the fact that the dark matter density field arising from 
Gaussian initial conditions exhibits a finite range of den¬ 
sities which occur with a significant probability. Note that 
this lower bound indeed agrees with the points obtained 
from numerical simulations shown in Fig.[t| On the other 
hand, the high-density bound p c is due to gravitational 
shock-heating which becomes the dominant energy source 
for dense regions. 

Secondly, there is a “warm” IGM phase (T ~ 10 4 —10 7 
K) which describes the gas heated by shocks arising from 
the gravitational energy of just collapsing objects. Because 
of the stochastic character of this energy source there is a 
broad scatter for this component around the “Equation of 
State” we have derived. This gas is restricted to a specific 
allowed region in the (p, T) phase-diagram. This expresses 
cooling and heating constraints as well as the properties of 
the underlying dark matter density held. We nevertheless 
obtained a curve, Eq.(p3|), which follows the mean trend 
of this “warm” phase. Its low density bound is set by the 
transition near p c with the “cool” IGM phase dominated 
by radiative heating from the UV background. The high- 
density bound is given by the intersection with the cool¬ 
ing curve (where bremsstrahlung cooling becomes domi¬ 
nant) or the density threshold A v ; r (beyond this point we 
have groups or clusters of galaxies). Note that these results 


are consistent with the calculations of Nath fe Silk (200T)| 
based on the Zel’dovich approximation. 

We can note that our results shown in Fig. [I] agree 
reasonably well with the outcome of numerical sim¬ 
ulations as displayed in Fig.11 in Dave et al. (19991 
(also shown by the points in our Figjlj) and Fig.9 in 
[Springel fe Hernquist (2001) . There is a small offset at 
low redshift for the normalization of the equation of state 
( pT] ) of the “cool” IGM and for the cooling region de¬ 
fined in eq.([|). Note that the latter could be remedied 
by adjusting the ratio t c ,ef/ti? which we simply set equal 
to unity in eq.(Q). Similarly, we could obtain a better fit 
to the numerical results for the “warm” IGM by tuning 
the r.h.s. in eq.(p5|) which defines the mass Mq. However, 
our goal is not to get the best fit to a specific numerical 
simulation (which would be of little value) but to explain 
the physics of the IGM. Moreover, Fig.^j shows that the 
simple procedure detailed in the previous sections already 
provides a good qualitative and quantitative description 
which should be sufficient for most purposes. Besides, as 
explained above it should be quite robust. In particular, it 
could be readily used with any cosmological parameters. 


2.6. Redshift evolution 


Finally, the four diagrams displayed in Fig .|l] show the evo¬ 
lution with redshift of the properties of the IGM, from 
z = 3 down to z = 0. We can see that the pattern does 
not evolve much qualitatively although the curves exhibit 
a quantitative shift with z. The characteristic tempera¬ 
ture of the “cool” IGM decreases slightly with time, in 
agreement with eq.(|2l|) since the density declines faster 
than 1/fir. Indeed, the cooling rate due to the expan¬ 
sion of the universe scales as 1/trr while the recombina¬ 
tion rate (which yields the density of neutral hydrogen in¬ 
volved in radiative heating) scales as p. On the other hand, 


the cooling constraint defined by t ce { = f# in Sect.2.2.1 
does not evolve much with time. In particular, its lower 
branch at T ~ 10 4 K is set by the atomic physics of hy¬ 
drogen ionization. Next, the characteristic temperature of 
the curve p+ t \ oc (T) associated with shock-heating due to 
gravitational clustering grows with time. This expresses 
the fact that the virial temperature associated with larger 
scales which turn non-linear later is higher. For instance, 
the velocity dispersion associated with galaxies is of order 
~ 200 km/s while for clusters it is ~ 1000 km/s. Following 
this evolution, the mean curve T g h(p) which describes the 
“warm” IGM enters the cooling region at high z while at 
low z it first crosses the high-density threshold A v ; r . This 
actually expresses the fact that at high z the collapsed 
halos built by gravitational clustering form galaxies since 
the gas undergoes a very efficient cooling, while at low z 
typical just-virialized halos are clusters which remain hot 
over a Hubble time and are still strong X-ray emitters 


(e.g., Valageas & Schaeffer (2000)). 
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Fig. 2. The phase-diagram of cosmological baryons at z = 3 and z = 0. As in Fig.^Jthe curves at <5 < A v ; r (z) are the 
equations of state of the “cool” and “warm” phases of the IGM. The dashed curve with a vertical part at A v j r (z) ~ 200 
and a branch towards higher 5 and T corresponds to galaxies. The two upper parts of the vertical line at A v j r (z) are 
groups (solid line at z = 0) and rare clusters (dotted line). The upper right dot-dashed line corresponds to the cool 
cores of groups and clusters where cooling has had time to develop. 


2.7. The IGM versus virialized halos 

As explained in the previous sections, we restrict the 
“mean equation of state” of the “warm” IGM to density 
contrasts below the threshold A v j r . Indeed, we consider 
that larger densities correspond to galaxies or clusters, 
that is virialized halos which are not part of the diffuse 
IGM. We display in Fig.|] the location in the (p, T) phase- 
diagram of these collapsed objects and of the IGM at red- 
shifts z = 3 and z = 0. 

Let us first consider the right panel, obtained for z = 0. 
The dashed line at S ~ 0 is the “cool” IGM while the 
solid curve at S c < 5 < A v j r is the “warm” IGM. The 
equations of state of these two components were derived 
in the previous sections and these curves are identical to 
those shown in Fig.[j]. For clarity we do not plot in Fig.|] 
the boundaries discussed in Fig. [l] which constrain the large 
scatter of the “warm” IGM phase. 

Next, the vertical line at <5 = A v j r (z) ~ 200 shows 
the overall density contrast associated with just-virialized 
objects. We divide this line into three parts. 

The first is in the low temperatures T ~ 10 5 K region. 
These objects with S = A v j r are located within the cool¬ 
ing region shown in Fig||. Hence their cooling time f c>e f 
is smaller than the Hubble time tn, see eq.([|). Therefore, 
within these halos the gas undergoes a very efficient cool¬ 
ing which leads to the formation of stars. As a conse¬ 
quence, these objects are small galaxies which are just 
being formed. This part of the vertical line S = A v ; r which 
is enclosed within the cooling region is shown by the lower 
vertical dashed line in Figj^. 

The remaining high temperature region may still be 
subdivided. The halos at S = A V i r which have a higher 
temperature do not cool over a Hubble time since they 
obey t Cj6 f > tjj (except in their center). Hence they cor¬ 
respond to groups or clusters which still contain hot gas 


which can be observed through its X-ray emission. One 
may divide this part into two components: groups ver¬ 
sus clusters. We identify groups with the low tempera¬ 
ture halos located below the cutoff p+,i oc shown in Fig.0. 
Therefore, they correspond to the typical just-virialized 
halos which form at z = 0. On the other hand, we identify 
clusters with the high temperature halos located above 
the cutoff p+.ioc- Hence they are rare massive objects 
which probe the high-density tail of the pdf V(pr). In 
other words, they correspond to the high-mass tail of the 
mass function. Thus, our distinction between groups and 
clusters is only based on the abundance of these objects 
(whether they correspond to rare or typical density fluc¬ 
tuations). However, we shall explain in a future paper 
(Valageas et ah (2002)) that this subdivision also marks 
the boundary between objects which, depending on the 
depth of their gravitational potential, are affected or not 
by a preheating of the gas through the gravitational pro¬ 
cesses which yield the IGM. Groups (resp. clusters) are 
shown in Fig.| by the vertical solid (resp. dotted) line. 
Thus, these three parts of the vertical line at 6 = A V i r de¬ 
scribe the location of small galaxies, groups and clusters. 

For completeness, we must point out that the lower 
part of the line at 8 = A v i r with T <, 2 x 10 5 K does 
not give all galaxies. Indeed, it is clear that there are 
some galaxies with a higher temperature. As discussed in 
Valageas fe Schaeffer (1999) these objects must still sat¬ 
isfy a somewhat stronger cooling constraint: 


^c,ef — tr[{z form) 


(26) 


where Zf or m is the redshift of formation of the object. Halos 
with S = A v ; r (z) simply have Zf or m = z (i.e. they are 
just being virialized) while higher density contrasts corre¬ 
spond to older objects which collapsed at earlier times. 
The constraint (|2^) ensures that cooling was very effi¬ 
cient when these halos formed so that their gas was able 
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to condense fast enough in order to be left unperturbed 
by the next generation of forming structures. This al¬ 
lows these objects to retain their individuality, stars and 
possibly a gaseous disk, even if they become embedded 
within larger structures (e.g., filaments or clusters) at later 
times. The condition (|2(j) at <5 > A v i r is shown by the 
dashed line in Fi g.|2j, where we use d the approximation 
tfiiz form) = t h \/(1 T A v i r )/(1 + d) (i.e. form) scales 
as the dynamical time of the halo). Thus, we are led to 
consider that the properties of these galactic halos (mass 
and size) do not evolve much after they form and reach the 
virialization density contrast 5 = A V i r . This method yields 
predictions for the properties of galaxies (luminosity func¬ 
tion, masses, metallicity,...) which agree with observations. 
Therefore, galaxies are defined by the dashed line which 
shows a first vertical part at 8 = A V j r and T ^ 2 x 10 5 
K (small galaxies which are being formed) and a second 
part at higher densities and temperatures (older and more 
massive galaxies). We can check that we obtain at T = 10 6 
K a radius R ~ 50 kpc which is indeed the typical size of 
present galaxies. By contrast, a halo with the same virial 
temperature at 6 = A v j r would give a size ~ 300 kpc 
(and a mass ~ 1O 13 M 0 ) which is uncomfortably large for 
a galaxy. This actually corresponds to a small group of 
galaxies. 

Finally, we note that the upper part of the cooling 
curve (0) shown in Fig.jl] at 6 > A v ; r describes the cool 
cores of groups and clusters. Indeed, it corresponds to 
high-density regions which have had time to cool today. 
This yields the core radius of present groups and clusters 
where cooling has just had time to come into play and 
possibly induce cooling flows. For instance, at T = 10 7 K 
we obtain a radius i? coo i ~ 150 kpc, which agrees with 
observations, while the virial radius is R V1I ~ 1 Mpc. We 
display this curve as the dot-dashed line in Fig.0. Note 
that this assumes that baryons follow the dark matter. 
However, at low z this is not necessarily the case since 
some “preheating” may modify the physics of the gas 
in low-temperature clusters (i.e. groups), see for instance 


Valageas fe Silk (1999)| and |Valageas fe Schaeffer (2000) 


This is not important for our purpose here which is mainly 
to describe the physics of the IGM. A specific study de¬ 
voted to this problem and the entropy of the gas in the 
light of the results described in this paper is presented in 


Valageas et al. (2002) 


Note that for these collapsed objects (galaxies, groups, 
clusters and cool cores) the density contrast 6 and the 
temperature T shown in Fig.0 only refer to the mean 
density contrast and virial temperature of the halo over 
the relevant radius R (which may be different from the 
virial radius for galaxies and cool cores). At the center of 
the halo the gas density is larger. Moreover, for galaxies 
and cool cores the gas temperature can be significantly 
smaller since cooling is very efficient (and we should 
take into account feedback from supernovae). Therefore, 
contrary to the IGM where the location in the (p, T) 
phase-diagram directly gave the properties of the gas, 
here the ( p,T ) plane only shows the overall properties 


of the dark matter halos associated with each object. 
Finally, we must point out that these classes of objects 
are not exclusive of each other. Indeed, cooled cores are 
obviously embedded within groups and clusters while 
galaxies can be found within filaments, groups and clus¬ 
ters. Therefore, some of the mass associated with galaxies 
or cool cores is counted within the matter attached 
to filaments, groups or clusters (this can be handled 
using the methods of Valageas & Schaeffer (1997) 


see 


and 


Valageas & Schaeffer (1999) for galaxies 
|Valageas fe Schaeffer (2000)| for clusters). 

Thus, as explained above the (p, T) diagram shown in 
Fig.0 describes the physics of cosmological baryons and 
their distribution between different phases. 

The redshift evolution up to z = 3 can be easily de¬ 
rived from Fig.0 and for illustration we show in the left 
panel in Fig.0 our result at z = 3. The location of the vari¬ 
ous curves evolves as explained in Sect. |2.6| but we can also 
note some qualitative changes. Firstly, with our definitions 
we see that there are no more groups. This means that 
the typical objects which collapse at z = 3 form galax¬ 
ies as they exhibit efficient cooling. To be more precise, 
one can still find clusters along 6 = A v j r , at high temper¬ 
atures where cooling is not very efficient. However, they 
are located within the far tail of the mass function and 
they correspond to extremely rare events. Secondly, we 
note that the curve which describes the mean “equation 
of state” of the “warm” IGM stops below A v ; r . Indeed, 
as seen in Fig-0 its high-density end-point is now given 
by the intersection with the cooling region. This is merely 
another consequence of the fact that we typically form 
galaxies and not groups (this feature coincides with the 
“disappearance” of groups). At an even higher redshift 
z ~ 4 the “warm” IGM component almost disappears as 
the equation of state of the “cool” IGM extends up to 
the cooling region. However there still remains a “warm” 
phase because the allowed region in the (p,T) plane has 
not vanishedj 2 |. At these high redshifts most of the mass 
is within a roughly uniform “cool” phase. We still have 
some rare collapsed halos which correspond to galaxies 
and some non-linear structures where gravitationally in¬ 
duced shocks heat the gas but these latter regions are 
severely restricted by the high efficiency of cooling pro¬ 


cesses. 


3. Distribution of matter 

3.1. Mass fraction within different phases 

Thus, in Sect.0 we have discussed how the properties of 
cosmological baryons could be seen through a (p, T) phase- 
diagram. In particular, we distinguished a “cool” and a 
“warm” IGM phase. In this section, we derive the fraction 
of matter enclosed within these various components. 

2 These objects are not yet virialized and still in their col¬ 
lapse phase: their gas is still expected to undergo shocks which 
lead to the “warm” IGM, but the latter then rapidly cools, so 
its physics is rather different. 
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Let us first discuss the “cool” IGM phase. As described 
in Sect j2..'i| this corresponds to the Lyman-a forest, that is 
moderate density fluctuations governed by the ionization 
and the heating due to the UV background. Since its tem¬ 
perature is not zero this gas probes the dark matter den¬ 
sity field over the scale R defined in eq. (||), which describes 
the length-scale over which pressure can homogenize the 
baryonic matter distribution. Note that this length de¬ 
pends on the temperature since R oc y/T. Then, we wish 
to express the fraction of matter within the “cool” IGM in 
terms of the pdf V(pr) over the scales R associated with 
these clouds. To do so, we first note that these scales are 
within the non-linear regime as shown by the curve T n i(p) 
in Fig.[l] which marks the transition to the linear regime. 


Then, as discussed in [Valageas fe Schaeffer (1997)| , in the 


highly non-linear regime the pdf V{pr) shows the scaling: 


V(pr) = =2 M x ) with x = 


(27) 


which defines the scaling function h(x). This holds for den¬ 
sities larger than the low-density cutoffs p_ discussed in 
Sect. 2.2.4. Note that if the stable-clustering ansatz is valid 
this function h(x) does not depend on redshift, as long as 
one remains in the highly non-linear regime. Then, the 
fraction of matter enclosed within spherical cells of radius 
R with a density between the thresholds pi and p 2 is: 


F m = 


rP 2 r x 2 ( \ r 

/ d pr P rV(pr)= / — x 2 h(x). 

J pi Jx 1 ^ 


(28) 


The advantage of the scaling variable x is that the 
last term in eq.(p8|) still provides an estimate of the 
fraction of matter F m associated with a curve (p, R) 
wh ich is not necessarily at con stant radius, as discussed 
in Valageas fe Schaeffer (1997) T his prediction is com- 
pared with numerical simulations in Valageas et al. (2000) 
where we investigate the two cases of a constant density 
and a constant radius. Therefore, in order to estimate the 
fraction of matter associated with the Lyman-a forest we 
use the variable x and we write: 


F m — J -j- x 2 h{x ; i?, z), 


with 


x h(x\ R, z) = PrP{pr ; R, z) 


(29) 


(30) 


where V{pr; R, z ) is the value of the pdf V(pr) at scale R 
and redshift z and the integral is taken along the relevant 
curve ( p,R ). The pdf V{pr) is obtained from eq.(A.2) 
where the parameter n which enters eq.(A.5) is computed 
at the relevant scale R. Thus, the expression ( |29|) w ould 
be exact within the framework described in App jA.lj if the 
Lyman-a forest were defined by a constant scale R. The 
use of the variable x in eq. (|29|) is only meant to handle 
the fact that R actually varies with p. In particular, we 
must point out that our procedure ([29|) with the definition 
( lso| ) does not assume that the stable-clustering ansatz is 
valid nor that we consider highly non-linear scales. Indeed, 



Z 


Fig. 3. The distribution of baryonic matter. The dashed 
line which increases with redshift is the fraction of matter 
within the “cool” IGM phase, as computed from eq.(p9|) 
and eq.(pl|). The dot-dashed line gives the mass fraction 
associated with the “warm” IGM component. The solid 
line is the fraction of matter within collapsed objects. The 
filled symbols show the results of the numerical simula¬ 
tions from Dave et al. (2001) (panel D2 in Fig.l) for the 
“cool” IGM (squares), the “warm” IGM (triangles) and 
condensed gas (circles). The empty squares show the re- 
sults of the simulation presented in Dave et al. (19991 for 
the “cool” IGM phase (Fig.12). 


the approximation (^9|) merely expresses the mass fraction 
associated with certain clouds of size R and density pr in 
terms of the pdf V(pr', R, z) measured at this point p#, R). 

As described in App jA.l] our approximation for the 
pdf applies both to the linear and non-linear regimes. 
Moreover, this formalism does not imply that these clouds 
are virialized objects which have reached an equilib¬ 
rium state. The “extended” function h(x: R, z ) defined in 
eq.(|3^) depends on scale and redshift. It is merely an¬ 
other way to write the pdf V(pr; R, z) (which goes over 
to the one relevant for the highly non-linear regime in this 
limit and is thus used by analogy). An important point of 
eq.(||^) is that it yields the mass fraction associated with 
various objects from the pdf V(pr) of the non-linear den- 
sit y field. This is a strong adv antage since, as described 
in Valageas & Schaeffer (1997), one can show that in the 
highly non-linear regime, assuming the stable-clustering 
ansatz is valid, the cloud-in-cloud problem can be han¬ 
dled in a satisfactory way. We can expect this property 
to extend to the regime we consider here. Note also that 
eq.(|2^) allows us to count objects which are not necessar¬ 
ily defined by a constant density contrast. 

Thus, the integration of eq. (E9) along the Equation of 


state (|21| ) yields the fraction of matter Fm,cooi within the 
“cool” IGM phase. Next, we could try to use the same 
method for the “warm” IGM, that is to integrate eq.(^9|) 
along the mean “equation of state” (^). However, this 
procedure is not satisfactory since the filaments and large 
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non-linear structures also contain galaxies and groups 
which should be subtracted. As a consequence we rather 
compute the mass fraction associated with virialized halos 
at S = A v ; r . That is, we integrate eq.^) along the ver¬ 
tical line at A v j r shown in Fig| starting from the lower 
point which is given by the lower intersection with the 
cooling curve. This yields the mass fraction F m , v - lr associ¬ 
ated with virialized halos (galaxies, groups and clusters). 
Then, we simply define the mass fraction A rn . warrn within 
the “warm” IGM by: 


F, 


m, warm 


= l-F 


m,cool 


— Fry 


(31) 


so that the sum is unity. This provides the distribution of 
matter within these three components: the “cool” IGM, 
the “warm” IGM and collapsed halos. 

We show our results in FigJU The Lyman-cc forest cor¬ 
responds to the dashed line which grows at larger red- 
shift. Indeed, at higher z gravitational clustering was less 
advanced so that a smaller fraction of matter had been 
shock-heated to high temperatures or embedded within 
collapsed objects through the building of large scale struc¬ 
tures. Note that the fractions of matter within collapsed 
halos and the “warm” IGM are of the same order. Indeed, 
both components are related to the formation of non¬ 
linear gravitational structures. As noticed in Sect. 2.7 we 
obtain Am,warm = 0 at 2 > 4. This actually means that 
the mass fraction within the “warm” IGM is very small 
and within the inaccuracy of our computation. 

We can check that our results agree with the 
outcome of the numerical simulations described in 


Dave et al. (1999) and Dave et al. (2001) . The mass frac¬ 
tions we obtain for the “warm” IGM and virialized ob¬ 
jects are also in reasonable agreement with simulations 
although it is difficult to make a detailed comparison. 
Indeed, our separation between different components is 
not exactly the same as in the simulations^ and the pre¬ 
dictions of the various simulations exhibit a significant 
scatter as shown by the comparison between the filled and 
empty squares in Figg (see also the various panels in Fig.l 
in Dave et al. (2001 )| ). 


3.2. Temperature distribution of the warm phase 

We have shown in Figg the fraction of matter within the 
“warm” IGM. It is also of interest to consider more pre¬ 
cisely the temperature distri but ion within this “warm” 
IGM. As we explained in Sect.|d], the mass function which 
describes the “warm” IGM cannot be obtained by directly 
integrating eq. (f29|) along the mean “equation of state” 

3 We define, in what we consider a somewhat more logical 
way, the various components according to their equation of 
state in the phase-diagram, which is easily recognizable for 
the “cool” or “warm” IGM, or for the virialized halos. This 
definition, as expected, changes accordingly with redshift. In 
the simulations, on the other hand, a simple temperature cut¬ 
off is used. For the components containing a small fraction of 
the baryons, this may introduce some difference. 



4 5 6 7 e 


log(T) [K] 


Fig. 4. The temperature distribution p(T) of the “warm” 
IGM component at redshifts z = 2 and z = 0 (normalized 
to unity). The mean temperature grows with time as larger 
scale structures form. 


because the regions which are about to collapse also 
contain smaller virialized halos (galaxies, groups). Thus, 
we first need to subtract these substructures which yields 
eq.(l3l|) for the mass fraction. This subtraction is more 
difficult when we consider the detailed temperature dis¬ 
tribution within the “warm” IGM but we shall assume 
that it can be described by a simple overall renormaliza¬ 
tion. This would be exact if the fraction of matter enclosed 
within such virialized objects would be the same for all 
“warm” IGM regions (i.e. filaments). In fact, we can ex¬ 
pect some correlation between these substructures and the 
larger-scale filaments but this procedure can be seen as a 
zeroth-order approximation. As a consequence, following 
eq.(|28|) we write the fraction of matter p(T)dT/T within 
the range [T, T + dT] as: 


« x 2 h(x) —, 

1 X 


(32) 


which we compute along the curve ( f23| ) which describes 
the mean “equation of state” of the “warm” IGM. In 
eq.(|3^) and in the following we do not write explicitly 
the time and scale dependence of h(x\R,z). The nor¬ 
malization constant in eq.([g|) is set by the constraint 
f p(T)dT/T = F m , warm , which is the total mass fraction 


within the “warm” IGM derived in Sect.3.1. 

We display our result normalized to unity in Figg 
Of course, we recover the fact that the mean tempera¬ 
ture increases with time as larger non-linear structures 
build up. This could also be seen from Figg. In partic¬ 
ular, the mean trend curve of the “warm” IGM in Fig jg 
at z = 0 is very steep because the characteristic temper¬ 
ature of the non-linear structures which collapse at that 
time is of order T co \ ~ 2 x 10' K which is much larger 
than the temperature reached by radiative heating from 
the UV background T a ~ 5 x 10 3 K. Therefore, gravita¬ 
tional shock-heating easily heats the gas up to T <; 10 6 K. 
However, note that there is still a significant fraction of 
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the matter at lower temperatures T < 10 6 K, in agreement 
with Dave et al. (2001) . Although the procedure (Js^) is 
only a simple approximation, the comparison of Fig.ETwith 
Fig.5 in Dave et al. (2001) shows that it captures the main 
trend. Of course, as explained in Sect.2.4 the “warm” IGM 
shows a broad scatter in the (p, T) plane so that the ac¬ 
tual boundaries of the temperature distribution are not 
as sharp as in Fig.|j. On the other hand, note that in our 
description we also have some matter at T > 10' K but 
we associate this gas with massive X-ray clusters. In any 
case, our result displayed in Fig^ appears to provide a 
reasonable description of the “warm” IGM properties. 


4. Warm IGM two-point correlation 

Another check on our description of the “warm” IGM com¬ 
ponent is provided by its two-point correlation function 
£ w (r) which measures its spatial clustering. Hence in this 
section we evaluate £ w (r), which is defined by: 


1 + £w(r) 


(p w (ri)p w (ri + r)) 

<Pw > 2 


(33) 


where p w (r) is the density of the “warm” IGM component 
at point r. In eq. (|33|) the averages (..) are over the realiza¬ 
tions of the density fields or over space (assuming ergod- 
icity). First, following the procedure described in Seet. [S.2| 
we can write the mean “warm” IGM density as: 

/ d(2? ~o 

- x2 h(x)-T7Vp = P-Fm.warm (34) 

x M 

where M, R and V are the mass, the radius and the vol¬ 
ume associated with the point (p, T) in the phase-diagram. 
The normalization constant A f is defined by: 


/ dx 

x h{F) = F m warm . (35) 

x 

Note that we always have (p w ) < p since the “warm” IGM 
component does not contain all the baryonic matter. Next, 
the second-order moment is given by: 


/ dx 

—-x 2 ft.(xi)p w (r, x\) 

Xi 


(36) 


where p w (r, xi) is the mean “warm” IGM density at 
distance r from a point embedded within the “warm” 
IGM component with the parameter X\. Then, as in 
Valageas et al. (2001a)| , we write p w (r, Xi) as the sum over 
two contributions. First, there is a probability F{r/R\) 
that the second point at iq -|- r belongs to the same region 
of thickness Ri of the “warm” IGM. Assuming spherical 
regions we have: 


S )j> p 3 

r<R 1 :F = l- — + ^,r>R 1 :F = 0. (37) 

This gives the first contribution pw^(r, Xi) as pi) 1 (r, x±) = 
p\F(r/R\). Second, we write the contribution p^ (r, xi) 



log(r) [Mpc] 


Fig. 5. The “warm” IGM two-point correlation £ w (r) at 
redshift z = 0 (dashed line). The solid line shows the two- 
point correlation £(r) of the dark matter density field. 


to p w (r, xi) of other regions of the “warm” IGM compo¬ 
nent with parameter X 2 as: 


Pw\ r ,xi) = [1 - F{r/Ri)} Mp /— x 2 ft,(x 2 ) 

J x 2 

x [1+ 6(xi)&(x 2 )£(r)]. (38) 


Here £(r) is the two-point correlation of the underlying 
dark matter density field while 6(xi)6(x 2 ) is the bias as¬ 
sociated with the regions x\ and x 2 . This bias is com¬ 
puted from the generating function introduced in eq.(A.2) 
as described in Bernardeau fe Schaeffer (1992) (see also 
jValageas et al. (2001b) for a comparison with observa¬ 
tions). Thus, we write the second-order moment as: 


(Pw(ri)Pw(ri + r)) = Mp 2 J ^-x\h{xi)^p R1 F(r / Ri) 

+ (1 - F(r/Ri)) M J ^p-x 2 /i(x 2 ) [1 + &(xi)6(x 2 )£(r)] j 

(39) 


This yields the “warm” IGM two-point correlation £ w (r) 
through eq. (|33|) . 

We show our result at z = 0 in Fig.||. On large scales 
the “warm” IGM correlation closely follows the dark mat¬ 
ter correlation because i) we assumed that baryonic mat¬ 
ter follows the dark matter density field smoothed over 
the scale R which describes the processes at work (R is 
set by the location on the (p, T) plane) and ii) the “warm” 
IGM component is mainly made of typical density fluctu¬ 
ations. This is seen by the fact that, by definition, most of 
its matter is within the allowed region in the (p, T) phase- 
diagram, below the high-density cutoff p+ : \ oc (T). This im¬ 
plies that the biases which enter the expression ( |39| ) are 
close to unity since we do not select rare high densities 
or voids by looking at the “warm” IGM phase. In fact, as 
seen from Fig.^j we obtain a mean bias which is slightly 
smaller than unity. However, this does not take into ac- 
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count the broad scatter in the (p, T) plane of the “warm” 
IGM component. 

On the other hand, at small scales r ^ 300 kpc the two- 
point correlation function £ w (r) flattens and it reaches a 
finite limit at r = 0. This expresses the fact that the non¬ 
zero temperature of the gas reached by radiative heat¬ 
ing and shock-heating (T > 10 4 K) gives rise to pressure 
effects which homogenize the baryonic distribution over 
some scale R. A lower limit to this length R is set by the 
scale R c associated with the critical point (p c ,T c ) in the 
phase-diagram. Indeed, as shown by the curved solid line 
in Fig-1 at redshift z = 0 the length scales associated 
with the “warm” IGM run from R c up to scales of order 
Ro which are turning non-linear. This sets the knee of the 
two-point correlation £ w (r) at r ~ 300 kpc at z = 0. As no¬ 
ticed in Sect.2d, the result displayed in Fig.|| shows that 
the characteristic thickness of the filaments which build 
the “warm” IGM is of order 300 kpc at z = 0. We find 
this, at a given redshift, to be indeed much smaller than 
the scale which is just turning non-linear, in agreement 
with numerical simulations. 

The value at r = 0 yields the clumping factor C w which 
we define by: 


C w = 


(Pi) 

(Pw > 2 


= 1 +£w(0). 


(40) 


We obtain C w ~ 160 at z = 0. Finally, we can check that 
our results agree wit h the outcome of th e numerical sim¬ 
ulations described in Dave et al. (2001) (see their Fig.7). 
In particular, we recover the knee at r ~ 300 kpc and a 
reasonable value for the clumping factor C w . Note how¬ 
ever that the dispersion in the (p, T) phase-diagram may 
actually lead to a slightly higher value. 


5. Soft X-ray background 

Finally, we check in this section whether the soft X-ray 
background emitted by the “warm” IGM in our model is 
consistent with observations. Experiments show that the 
extragalactic soft X-ray background flux in the 0.1 — 0.4 


keV band is of order 20 — 35 keV cm 


-2 


sr 


keV 


-l 


(e.g., |Warwick fe Roberts (1998)|) . However, most of this 
flux is due to AGN so that the contribution from diffuse 
gas should obey the constraint F v <, 4 keV cm -2 s _1 sr -1 
ke V -1 in this frequ ency range (e.g., see the discussion 
in Wu et al. (2001)). Using the Press-Schechter formal¬ 
ism ( Press fe Schechter (1974) ), assuming that the X-ray 
emitting gas is embedded within just-virialized halos with 
a density contrast A c {z) ^ 177, some previous studies ob¬ 
tained a flux which is higher than this upper bound (e.g., 
Pen (1999), Wu et al. (2001)). This led these authors to 
infer that a non-gravitational heating source is needed to 
unbind the gas from groups and small clusters in order to 
reduce their contribution to the soft X-ray background. 
However, the diffuse “warm” IGM component is made of 
filamentary structures with density contrasts lower than 
A c (z). This gives a smaller X-ray background since the X- 
ray emission is proportional to the squared density of the 


gas. Moreover, the fraction of matter within this “warm” 
diffuse phase is sm aller than 30%, see Fig|| As argued in 
Dave et al. (2001)| this reduced X-ray background could 
then be consistent with observational constraints. This 
conclusion also agrees with the outcome of numerical sim¬ 
ulations (Croft et al. (2001), Phillips et al. (2001 )| ). 

Therefore, we compute here the soft X-ray background 
due to the diffuse “warm” IGM phase, in the 0.1 — 0.4 
keV band. Note that the X-ray background due to re¬ 
solved sources ( AGN, clusters, cooling gala xies) was al¬ 
ready studied in Valageas fc Schaeffer (2000) from an ana¬ 
lytic model similar to the one used here to describe the un¬ 
derlying dark matter density field. The X-ray flux F Vl _„ 2 
within this frequency band can be written: 


F v ^ vn = 7.13 x 10" 34 keV cm 


2 ° 1 sr 1 


x [ dz Al + z) 2 Af(z) [ — x 2 h{x)(l+6)T 1 / 2 
J dz J x 

X (e~ hv 1 ( 1 + 2 )/( fe:r ) _ e -hv 2 (l +z)/(kT)'j , ( 41 ) 


where the temperature T is in Kelvin units. Here we 
used again eq.(|||) with the normalization Af(z) defined 
in eq.(^5[). The factor (1 + S) in eq.(ffd|) comes from the 
fact that the X-ray emission is proportional to the squared 
density n 2 . The integral over x in eq.(^) follows the 
curve T g h(/o), from eq.(^), in the (p, T) phase-diagram. 
Then, our calculation gives an X-ray background flux 
F Vl _„ 2 = 0.07 keV cm -2 s _1 sr -1 in the 0.1 — 0.4 keV 
band. This yields a mean differential flux F u = 0.25 keV 
cm -2 s _1 sr -1 keV -1 in this band. Thus, we can check 
that our result is consistent with observations. Note that 
the scatter of the “warm” IGM in the (p, T) plane should 
lead to a slightly larger flux. Therefore, this component 
might be observable through its X-ray emission in the near 
future. 


6. Conclusion 

In this article we have shown that the (p, T) phase- 
diagram of cosmological baryons can be understood in 
very simple terms. It is strongly model-independent since 
the relevant scales are set by atomic physics (hydrogen 
ionization) and the basic cosmological setting (the rate of 
expansion and the scale which marks the transition to the 
non-linear regime). Building a bridge between the mod¬ 
eling of galaxies, groups and clusters, as well as Lyman- 
a absorbers, we have constructed in this paper a phase- 
diagram which includes also the IGM, and thus provides 
a phase-diagram of all baryons. In this diagram, all these 
objects form well-defined populations. 

Firstly, we have distinguished a “cool” IGM phase cor¬ 
responding to the Lyman-a forest. It follows a well-defined 
equation of state on the (p, T) plane. We have shown that 
the properties of this gas do not depend much on its pre¬ 
vious history because of the form of the recombination 
coefficient. This explains why the scatter of this compo¬ 
nent obtained in numerical simulations is quite small. The 
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fraction of baryons in this components, as well as its evo¬ 
lution with redshift is seen to be consistent with earlier 
modelling of the frequency of occurrence of the Lyman-a 
lines. We however obtain (consistently with the numeri¬ 
cal simulations) 38% of the baryons in the cool IGM at 
z=0. This may be a little high with a mass fraction closer 
to 30% being consistent with observations. This must be 
considered as being within the error bars of our analytical 
calculation (and the error bars of the simulations!). 

Secondly, we find that a “warm” IGM phase is formed 
by the gas which has been shock-heated to larger temper¬ 
atures T ~ 10 4 — 10 7 K as non-linear gravitational struc¬ 
tures appear. The dependence on the stochastic gravita¬ 
tional potential entails a broad scatter in the (p, T) plane. 
However, we have explained that robust constraints only 
allow a closed region in the (p, T ) phase-diagram for this 
component. We have also defined a simple curve which 
represents its mean behaviour and plays the role of the 
“warm” IGM Equation of State. The latter is found to 
be consistent with the outcome of numerical simulations 
( Dave et al. (2001) ), but differs from the simple power-law 
fit of a mean curve done there. This is because its locus 
is based on physical considerations, which we have thor¬ 
oughly justified. The “warm” IGM is seen to be due to 
large structures on their way of collapsing under the ac¬ 
tion of gravity. The quite simple picture behind this model 
is that the collapse of an object turning non-linear (that 
is with average overdensity above 5) induces shocks which 
heat patches of much smaller size, to a temperature of the 
order of the locally available kinetic energy of the collapse, 
creating local conditions within hot spots close to hydro¬ 
static equilibrium. We are able to estimate this size, as 
well as the temperature of the latter, and to show they 
are in agreement with the simulations. Definitely, these 
objects are not just Lyman-a absorbers at larger densi¬ 
ties. The hot gas which lies within clusters and galactic 
halos, on the other hand, corresponds to the high-density 
“continuation” (e.g., S )> 200) of the “warm” IGM curve, 
in a region of phase-space where virialization insures that 
all the available gas has been gravitationally heated. 

Next, we find two quite different regimes for the 
“warm” IGM. At high redshift, some of the “warm” com¬ 
ponent enters the cooling region of the (p, T) plane: this 
gives rise to galactic disks and stars. Some of the baryonic 
matter, also, should lie on the low temperature branch of 
this cooling region, in agreement with the results of numer¬ 
ical simulations ( Dave et al. (1999)| ). At low redshift, the 
“warm” IGM no longer goes into the cooled region, and 
it provides the origin of the hot gas in clusters. The evo¬ 
lution with redshift of the haryon phase-diagram provides 
a natural confirmation that at high z collapsed halos form 
galaxies while at low redshift they build groups or clusters. 
The warm phase occupies 24% of the baryon fraction at 
z = 0 in our analytic model. 

Then, we have checked that our results for the fraction 
of matter enclosed within the various phases and the two- 
point correlation function of the “warm” IGM component 


agree with numerical simulations (e.g., Dave et al. (1999) 


and Dave et al. (2001)). This confirms the validity of our 
analysis. Note however that the latter relies on simple 
physical considerations and it is independent of the find¬ 
ings of the former. Besides, the soft X-ray background due 
to the “warm” IGM is consistent with the upper bound 
set by observations. Our prediction actually is not far from 
the observational limit: this offers the prospect of measur¬ 
ing this X-ray emission in future observations. 

We predict that approximately 60% of the baryons are 
accounted for at z=0 in the cool/warm phases of the in- 
tergalactic medium. The remaining baryons (i.e. 40%) are 
in collapsed structures. With a total of f \ ss 0.045, this 
amonts to flbcoii ~ 0.018. The observed luminous compo¬ 
nent of the baryons (stars, remnants and gas) is estimated 
to be about fl ga i ss 0.004, as noted in the introduction, and 
the hot gas in clusters represents fl g hot ~ 0.003. Therefore, 
there remains an “unidentified” or dark baryon compo¬ 
nent, flbdm ~ 0.011. This is at least twice the stellar- 
like component but only about 1/4 of the “identified” 
baryons (IGM, stellar components and very hot gas). In 
our model this dark baryonic matter corresponds to (cool 
?) gas within galactic halos and groups, which has not 
been observed yet. 


Acknowledgements. We thank R. Dave for prov iding us with 
some of the data published in Dave et al. (1999). 


Appendix A: Dark matter density field 

In this appendix we describe the analytic procedure we 
use to model the dark matter density field. Our method 
provides a simple means to estimate the pdf of the density 
contrast from the quasi-linear up to the highly non-linear 
regime, parameterized by only one variable (the skewness). 
This should be sufficient for our purpose which is merely 
to obtain the main properties of the IGM. Note that by 
construction the variance and the skewness of the pdf we 
obtain below agree with the results of numerical simula¬ 
tions, in both the quasi-linear and the highly non-linear 
regimes. 


A.l. The pdfV(pR) of the overdensity at scale R 

First, we note £ the moment of order two of the density 
contrast Sr over a spherical cell of radius R and volume V. 
It can be expressed in terms of the two-point correlation 
function £(xi,X2) = £(|xi — X2I) as: 


£=<<&> = 


IV 


dxx dx 2 
V 2 


£( xi , x 2 ). 


(A.l) 


The two-point correlation f(x) is the Fourier transform 
of the non-linear power-spectrum P(k ) which we compute 
from the linear power- spectrum P^(fc) using the analytic 
formulae obtained by Peacock fc Dodds (1996) from fits 
to N-body simulations. Of course, in the quasi-linear limit 
£ —> 0 we have £ = a 2 . Then, it is convenient to intro- 
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duce the generating function <p(y) defined by the inverse 
Laplace transform: 


/ +ioo 

-ioo 


d y 


— e 


[PRV-V>(.v)]/i 


27 Tit; 


(A.2) 


where V(pr) is the pdf of the overdensity over a cell of 
radius R: 


Pr = 1 + Sr. 


(A.3) 


= -y C(t) 


v{y) = y C{r) + 


(A.4) 


with: 


« t) = ( 1+ D 


This yields a family of functions ip(y) which depends only 
on the parameter k. 

Such a relationship arises naturally in the exact 
derivation of ip{y) in the qua si-linear limit (Bernardeau 
(1992,1994), Valageas (2002) ). Besides, it also appears 
in the non-linear regime within the framework of the 
“minimal tree-model” ( Bernardeau fc Schaeffer (1992) 


Schaeffer (1985) 

). It can be seen that the phenomenolog- 

ical form ( 
As in 

A.5 

) for C(t) also yields reasonable results. 

Bernardeau & Valageas (2000) 

, we choose the 


A.2. High and low density cutoffs p±(R) 

Now, we can estimate the high and low density cutoffs 
p±{R) of the pdf V( PR ) as follows. To each overdensity pr 
we asso ciate the Laplace variable y c {pR) (and r c through 
eq.(A.4)) defined as the saddle-point of the exponent in 
eq.(A.2). This can be written: 


p R = ip'(y c ) = C(t c ), 


(A.7) 


The advantage of the generating function <p(y) is that it 
reaches a finite limit in the quasi-linear regime £ <C 1. 
Moreover, in this regime one can derive <p(y) in a rigor¬ 
ous manner from the equations of motion (Bernardeau 
(1992,1994), |Valageas (2002) ). On the other hand, in the 
highly non-linear regime £ 1 the generating function 

(p(y) should also reach a finite (astonishingly close to the 
former, but seemingly not identical) limit if the stable¬ 
clustering ansatz is valid ( Balian fe Schaeffer (1989)] ) . 
This is seen to agree reasonably well with the results 
of numerical simulations (e.g., Valageas et al. (2000)|) . 
Therefore, hereafter we use the following simple param¬ 
eterization for the generating function ip(y): 


where we used eq.(A.4) to derive the last equality. Note 
that Sr and the pair (y Cl t c ) are of opposite signs. The 
reason for the introduction of y c is that for extreme over¬ 
densities pr —> oo or pr —> 0, that is in the tails of the pdf 
V(pr), the integral ( |A.2|) is dominated by the c ontribu- 
tion from y ~ y c (e.g., |Balian fe Schaeffer (1989) ). Then, 
it is clear that the density cutoffs of the pdf are given by 
the points p± where the argument of the exponential in 
eq. (0) at the saddle-point is of order unity and negative 
(note that it is zero at the mean overdensity (pr) = 1). 
Indeed, for larger <5/j| the argument becomes more neg¬ 
ative than —1 which leads to the exponential falloffs of 
V(pr). Therefore, we define the high and low density cut¬ 
offs p± by the condition: 


p±y± - v>{y±) 


= —l, 


(A.8) 


where y± = y c (p±) as given by eq.(A.7). This relation can 


be most conveniently written in terms of the variable r 
which yields: 


(A.5) T ± = =F 




(A.9) 


where we used eq.(A.5). In the highly non-linear regime 
eq. Q yields r + -C — 1. However, as discussed in 


Bernardeau & Schaeffer (1992) the implicit system (A.4) 
is singular at a finite value r s = —«/(2 + n) which leads to 
an exponential tail for the pdf V(pr) ~ e~ PR ^ Xs ^. Hence 
for —\/2£ < t s we take t+ = t s and p+ = x s ff Thus, we 
obtain for the high-density cutoff p+(R): 




so as to recover 


parameter n which enters eq.(H 
t he skewness S 3 = ( S%) c /l ; predicted by HEPT 

([Scoccimarro &; Frieman (1999)) in the non-linear regime. 


£<1: P+ = l + yj2£, £ > 1 : p+ = x s £ 

and for the low-density cutoff p_ (R): 


e « 1 : P- = 1 - , 


?» 1 : P- ~ C 


- re /2 


(A.10) 


(A.ll) 


More precisely, we use: 

- 34 - 314 — 2 n ) 

£ < 0.1 : S 3 = — — (n+ 3), ^ > 10 : S 3 = ^ + 2n+1 (A.6) 

and we interpolate between both values in the range 
0.1 < £ < 10. Here we note n the slope of the linear- 
power spectrum at the scale of interest. Note that th is 
procedure was seen in Bernardeau fe Valageas (2000) to 


match the results of numerical simulations with respect 
to weak gravitational lensing effects. Therefore, it should 
be quite sufficient for our present purpose which is mainly 
to investigate the main trends of the IGM properties. 


Note that in the quasi-linear regime we recover |<5±| ~ 
y/l = a as required for Gaussian initial conditions. Then, 
the density cutoffs p± are close to the mean: p± ~ (pr) = 
1 for £ <C 1. On the other hand, in the highly non¬ 
linear regime the cutoffs are displaced to very low den¬ 
sities (p_ <C 1: rare voids) and high densities (p+ 1: 

collapsed halos). 
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